Evaluating a potential technique with local optical flow vectors for automatic organ-at-risk (OAR) intrusion detection and avoidance during radiotherapy.

Various techniques of deep inspiration breath hold (DIBH) have been used to mitigate the likelihood and risk of exposing the heart, an organ-at-risk (OAR) for unintended radiation during left breast radiotherapy. However, issues of reproducibility of these techniques warrant further investigation into the feasibility of detecting the intrusion of an OAR into the treatment field during intra-fractional treatment delivery. The increase of high-dose, low-fraction radiotherapy treatments makes it important to immediately adapt treatment once an OAR is detected in the treatment field. This proof-of-concept implementation includes an algorithm that detects and tracks the motion at the edges of a treatment field and a control algorithm that adapts the treatment aperture according to the motion detected. In accordance to the AAPM Task-Group (TG-132) report, image registration techniques should be verified with virtual and physical phantoms prior to clinical application. Since most OARs move as a result of respiration-induced motion, we have used a lung phantom to generate images of a generic OAR intruding into a treatment field with known velocity. The phantom was programmed to move with sinusoidal and lung patient tumor motion patterns and the accuracy of intrusion tracking and MLC adaptation were benchmarked with the ground truth-programmed motion of the OAR. The motions were recorded with an electronic portal imaging device (EPID). An optimal cluster size of 9 × 9 motion vectors was found to provide the smallest average absolute position error of 0.3 mm. A strong linear correlation between the adapted MLC leaves and the actual OAR position was observed. The algorithm had a mean position tracking error of -0.4 ± 0.3 mm and a precision of 1.1 mm. It is possible to adapt MLC leaves based on the motion detected at the edges of the irradiated field, and it would be feasible to shield an unplanned intrusion of an OAR into the treatment field.

Reducing the tracking drift of an uncontoured tumor for a portal-image-based dynamically adapted conformal radiotherapy treatment.

Accurate tracking of organ motion during treatment is needed to improve the efficacy of radiation therapy. This work investigates the feasibility of tracking an uncontoured target using the motion detected within a moving treatment aperture. Tracking was achieved with a weighted optical flow algorithm, and three different techniques for updating the reference image were evaluated. The accuracy and susceptibility of each approach to the accumulation of position errors were verified using a 3D-printed tumor (mounted on an actuator) and a virtual treatment aperture. Tumor motion up to 15.8 mm (peak-to-peak) taken from the breathing patterns of seven lung cancer patients was acquired using an amorphous silicon portal imager at ~ 7.5 frames/s. The first approach (INI) used the initial image detected, as a fixed reference, to determine the target motion for each new incoming image, and performed the best with the smallest errors. This method was also the most robust against the accumulation of position errors. Mean absolute errors of 0.16, 0.32, and 0.38 mm were obtained for the three methods, respectively. Although the errors are comparable to other tracking methods, the proposed method does not require prior knowledge of the tumor shape and does not need a tumor template or contour for tracking. Graphical abstract.

Feasibility of predicting tumor motion using online data acquired during treatment and a generalized neural network optimized with offline patient tumor trajectories.

PURPOSE: The accurate prediction of intrafraction lung tumor motion is required to compensate for system latency in image-guided adaptive radiotherapy systems. The goal of this study was to identify an optimal prediction model that has a short learning period so that prediction and adaptation can commence soon after treatment begins, and requires minimal reoptimization for individual patients. Specifically, the feasibility of predicting tumor position using a combination of a generalized (i.e., averaged) neural network, optimized using historical patient data (i.e., tumor trajectories) obtained offline, coupled with the use of real-time online tumor positions (obtained during treatment delivery) was examined. METHODS: A 3-layer perceptron neural network was implemented to predict tumor motion for a prediction horizon of 650 ms. A backpropagation algorithm and batch gradient descent approach were used to train the model. Twenty-seven 1-min lung tumor motion samples (selected from a CyberKnife patient dataset) were sampled at a rate of 7.5 Hz (0.133 s) to emulate the frame rate of an electronic portal imaging device (EPID). A sliding temporal window was used to sample the data for learning. The sliding window length was set to be equivalent to the first breathing cycle detected from each trajectory. Performing a parametric sweep, an averaged error surface of mean square errors (MSE) was obtained from the prediction responses of seven trajectories used for the training of the model (Group 1). An optimal input data size and number of hidden neurons were selected to represent the generalized model. To evaluate the prediction performance of the generalized model on unseen data, twenty tumor traces (Group 2) that were not involved in the training of the model were used for the leave-one-out cross-validation purposes. RESULTS: An input data size of 35 samples (4.6 s) and 20 hidden neurons were selected for the generalized neural network. An average sliding window length of 28 data samples was used. The average initial learning period prior to the availability of the first predicted tumor position was 8.53 ± 1.03 s. Average mean absolute error (MAE) of 0.59 ± 0.13 mm and 0.56 ± 0.18 mm were obtained from Groups 1 and 2, respectively, giving an overall MAE of 0.57 ± 0.17 mm. Average root-mean-square-error (RMSE) of 0.67 ± 0.36 for all the traces (0.76 ± 0.34 mm, Group 1 and 0.63 ± 0.36 mm, Group 2), is comparable to previously published results. Prediction errors are mainly due to the irregular periodicities between cycles. Since the errors from Groups 1 and 2 are within the same range, it demonstrates that this model can generalize and predict on unseen data. CONCLUSIONS: This is a first attempt to use an averaged MSE error surface (obtained from the prediction of different patients' tumor trajectories) to determine the parameters of a generalized neural network. This network could be deployed as a plug-and-play predictor for tumor trajectory during treatment delivery, eliminating the need for optimizing individual networks with pretreatment patient data.